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Abstract 



We present measurements of the structure function F2 in e^p scattering at HERA 
in the range 3.5 GeV^ < Q2 < 5000 GeV^. A new reconstruction method has ahowed a 
significant improvement in the resolution of the kinematic variables and an extension of the 
kinematic region covered by the experiment. At < 35 GeV'^ the range in x now spans 
6.3- 10^^ < X < 0.08 providing overlap with measurements from fixed target experiments. 
At values of above 1000 GeV^ the X range extends to 0.5. Systematic errors below 
5% have been achieved for most of the kinematic region. The structure function rises as 
X decreases; the rise becomes more pronounced as increases. The behaviour of the 
structure function data is well described by next-to-leading order perturbative QCD as 
implemented in the DGLAP evolution equations. 
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1 Introduction 



Measurements of neutral current (NC) deep inelastic scattering (DIS) at HERA []T], have 
revealed a rapid rise of the proton structure function F2 as Bjorken-x decreases below 10~^. 
Extensions of these measurements to low have shown that this rise persists down to 
values as low as 1.5 GeV^ 0, |]. 

In this report we present a new measurement of F2 from a DIS event sample eight times larger 
than that used in our previous analysis The increase in statistics combined with a new 
method that provides a substantially more precise reconstruction of the event kinematics has 
allowed us to extend the range of the measurements in x by an order of magnitude and to 
decrease the systematic uncertainties of F2 by roughly a factor of three. 

New data on F2 with similar statistics in a similar kinematic range have recently been presented 
by HI i. 

2 Experimental conditions 

The data presented here were taken with the ZEUS detector at HERA in 1994. HERA operated 
with 153 colliding bunches of 820 GeV protons and 27.5 GeV positrons, with a time between 
bunches of 96 ns. Additional unpaired positron (15) and proton (17) bunches circulated, which 
are used to determine beam related backgrounds. The proton bunch length was approximately 
20 cm (r.m.s.) while the positron bunch length was negligible in comparison which, together 
with run-to-run variations of the mean interaction position, leads to a length of the interaction 
region of 12 cm (r.m.s.) centered around Z=+Q cm. Q The data of this analysis correspond 
to a luminosity of 2.50±0.04 pb"^. Approximately 5% of the proton current was contained 
in satellite bunches, which were shifted by 4.8 ns with respect to the primary bunch crossing 
time, resulting in a fraction of the ep interactions occurring at {Z) = +78 cm. 

A description of the ZEUS detector can be found in H, ^. The components used in this analysis 
are briefly discussed. The uranium-scintillator calorimeter (GAL) covers 99.7% of the total 
solid angle. It consists of the barrel calorimeter (BGAL) covering the range 36.7° < 6 < 129.1° 
in polar angle, the forward calorimeter (FGAL) covering 2.6° < 6 < 36.7° and the rear calorime- 
ter (RGAL) covering 129.1° <e < 176.2°. The FGAL and RGAL are divided in two halves to 
allow retraction during beam injection. Each calorimeter part is segmented in electromagnetic 
(EMG) and hadronic (HAG) sections. Each section is further subdivided into cells of typically 
5 X 20 cm^ (10 X 20 cm^ in the RGAL) for the EMG and 20 x 20 cm^ for the HAG sections. Each 
cell is viewed by two photomultipliers (PMTs). Under test beam conditions the calorimeter 
has an energy resolution of a/E = 18%/ ^ E{GeV) for electrons and a/E = 35%/ ^E{GeV) for 
hadrons. The timing resolution of a calorimeter cell is less than 1 ns for energy deposits greater 
than 4.5 GeV. In order to minimise the effects of noise due to the uranium radioactivity on 
the measurements all EMG(HAG) cells with an energy deposit of less than 60(110) MeV are 
discarded from the analysis. For cells with isolated energy deposits this cut was increased to 
100(150) MeV. 

""^The ZEUS coordinate system is defined as right handed with the Z axis pointing in the proton beam 
direction, and the X axis horizontal, pointing towards the centre of HERA. 
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The tracking system consists of a vertex detector (VXD) a central tracking chamber 
(CTD) 0, and a rear tracking detector (RTD) enclosed in a 1.43 T solenoidal magnetic 
field. The interaction vertex is measured with a typical resolution along (transverse to) the 
beam direction of 0.4 (0.1) cm. For high momentum tracks (p > 5 GeV) the extrapolated 
position on the inner face of the calorimeters is known with a typical resolution of 0.3 cm. 

The position of positrons scattered at small angles to the positron beam direction is measured 
using the small angle rear tracking detector (SRTD) which is attached to the front face of the 
RCAL. The SRTD consists of two planes of scintillator strips, 1 cm wide and 0.5 cm thick, 
arranged in orthogonal directions and read out via optical fibers and photo-multiplier tubes. 
It covers the region of 68 x 68 cm^ in X and Y and is positioned at Z = —148 cm. A hole of 
20 X 20 cm^ at the center of the RCAL and SRTD accommodates the beampipe. The SRTD 
signals resolve single minimum ionizing particles and provide a position resolution of 0.3 cm. 
The time resolution is less than 2 ns for a minimum ionizing particle. 

The luminosity is measured via the positron-proton bremsstrahlung process, ep — > ejp, using 
a lead-scintillator calorimeter (LUMI) |[T^ which accepts photons at angles < 0.5 mrad with 



respect to the positron beam direction. The LUMI photon calorimeter is also used to tag 
photons from initial state radiation in DIS events. It is positioned at Z = —107 m and has 



an energy resolution oi a / E = 18%/yi?(GeV) under test beam conditions. In its operating 
position, however, it is shielded from synchrotron radiation by a carbon-lead filter and has an 
energy resolution of ajE = 26.5%/ ^ E{GeV), as determined from bremsstrahlung data. The 
position resolution is 0.2 cm in X and Y. In addition, an electromagnetic calorimeter positioned 
at Z = —35 m is used for tagging positrons scattered at small angles. 



2.1 Triggering 

Events were filtered online by a three level trigger system At the first level, the events are 
selected by the logical OR of the following conditions (further details can be found in [[TTI ): 

• Total EMC energy deposit in the BCAL is greater than 4.8 GeV. 

• Total EMC energy deposit in the RCAL, excluding the region closest to the rear beampipe, 
is greater than 3.4 GeV. 

• In the RCAL, the isolated positron condition (ISO-e) is fulfilled. The ISO-e condition 
requires that the isolated EMC energy deposit be greater than 2.5 GeV and that the 
corresponding HAG energy be less than 0.95 GeV or no more than a third of the EMC 
energy. The above condition is ANDed with the requirement that the total energy deposit 
in RCAL EMC is greater than 3.8 GeV. 

Events where the positron is scattered in the direction of the FCAL are triggered efficiently, by 
the hadronic final state, with the above requirements. 

Backgrounds from protons interacting outside the detector were rejected using the time mea- 
surement of the energy deposits in upstream veto counters and the SRTD. The average trigger 
efficiency for events that pass the off-line selection cuts (see Sect. ^ is above 99%, as determined 
from independent triggers and MC simulation. 
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In the second level trigger (SLT), background was further reduced using the measured times of 
energy deposits and the summed energies from the calorimeter. The events were accepted if: 

6sLT = E^^(l - cos^.) > (24 - 2E^) (GeV) (1) 

i 

where Ei and 6i are the energies and polar angles (with respect toX = Y = Z = cm) of 
calorimeter cells, and is the energy deposit measured in the LUMI photon calorimeter. For 
perfect detector resolution, 6slt is twice the positron beam energy (55 GeV) for fully contained 
DIS events while for photoproduction events, where the scattered positron escapes down the 
beampipe, ^slt peaks at much lower values. Proton beam-gas events which originate from 
inside the detector have energy flows which are concentrated in the forward direction and so 
also have small values of Sslt- 

The full event information was available at the third level trigger (TLT). Tighter timing cuts 
as well as algorithms to remove beam-halo muons and cosmic muons were applied. The 
quantity 6tlt was determined in the same way as 6slt- The events were required to have 
Stlt > (25 — 2E^) (GeV). Finally, events were accepted if a scattered positron candidate of 
energy greater than 4 GeV was found. In total 900853 NC DIS candidates satisfied the above 
trigger conditions. 



3 Monte Carlo simulation 



Monte Carlo (MC) event simulation is used to correct for detector acceptance and smearing 
effects. The detector simulation is based on the GEANT program |]12[ and incorporates our 
understanding of the detector, the trigger and test beam results. Neutral current DIS events 
are simulated using the HERACLES program |13] which includes photon and Z° exchanges and 
first order electroweak radiative corrections. The hadronic final state is simulated using the 
colour-dipole model CDMBGF [|14| including all leading order QCD diagrams as implemented 



in ARIADNE M for the QCD cascade and JETSET |16] for the hadronisation. The ARIADNE 



model provides the best description of the observed DIS non-diffractive hadronic final state |T7 
Diffractive events which have been observed in the data 19 1 by the occurrence of a large 



rapidity gap in the detector are simulated within ARIADNE by assuming that the struck quark 
belongs to a colourless state having only a small fraction of the proton's momentum. The 
parameters of the model are adjusted to be consistent with recent ZEUS measurements 
The MRSA E^ parton density parameterisations, modified at low as described in 



|2 



are 



used. For the final acceptance corrections, the events are reweighted with the help of a NLO 
QCD fit to the data in an iterative procedure as described in Sect. ^. For systematic checks, 
a sample of events was also generated using the LEPTO Matrix Element and Parton Shower 
(MEPS) g final state simulation. 

The shape of the vertex distribution used in the simulation is taken from non-diffractive pho- 
toproduction events, since for these events the vertex reconstruction efficiency is found to be 
both high (~ 90%) and independent of the Z position of the interaction. 

The effects of the uranium radioactivity were simulated using distributions obtained from ran- 
domly triggered events. 
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A MC event sample corresponding to an integrated luminosity of 2.5 pb ^ for > 1.8 GeV^ 
was generated. This sample was supplemented by an additional sample equivalent to 2.5 pb~^ 
with g2 > 20 GeV^. 

The main source of background in the data comes from the few photoproduction interactions 
which lead to the detection of a fake scattered positron. 'Minimum bias' photoproduction 
events were simulated using PYTHIA with cross sections according to the ALLM parame- 
terisation p5|. Photoproduction events corresponding to an integrated luminosity of 250 nb~^ 



were generated with a photon-proton center-of-mass energy W > 190 GeV. Events with smaller 
W values do not contribute to the photoproduction background due to the requirement on 6 
(see Sect. 6). 



4 Energy and angle measurements 

In the determination of the DIS kinematics, the GAL energy deposits are separated into those 
associated with the identified scattered positron, and all other energy deposits. The latter is 
defined as energy from the hadronic system, or hadronic energy. The kinematics of the event 
is then determined from: 

• The energy {E'^) and polar angle (6^) of the scattered positron. 

• The hadronic energy in terms of 

6^^Y.iEh-pzh) (2) 

h 

and 



\l h h 

where the sums run over all GAL energy deposits not associated with the scattered 
positron. 

The hadronic energy flow is characterised by the angle 7^ defined by: 

COS7, = (4) 

In the naive quark parton model picture of DIS the angle 7^ corresponds to the polar angle of 
the struck quark. 



4.1 Positron identification and efficiency 

The positron identification algorithm (SINISTRA94) is based on a neural network using infor- 



mation from the GAL ||26[ . The network separates deposits in the calorimeter which are due to 
electromagnetically showering particles from those which are of hadronic origin, by assigning 
a 'positron-probability'. A cut on this probability allows a clean identification of the scattered 



4 



positron.0 The efficiency for finding the scattered positron was determined from MC simulations 
to increase from 80% at 10 GeV to greater than 99% for energies above 15 GeV. The efficiency 
was checked by using QED-Compton events and by comparing results from this algorithm with 
other positron finding algorithms [ij and found to be consistent with expectations. 

4.2 Positron angle measurement 

The scattering angle of the positron is determined either by measurement of its impact position 
on the GAL inner face together with the event vertex, or, for O^. < 135°, from the parameters 
of a reconstructed track matched with the positron impact position. 

The impact position of the scattered positron at the calorimeter is reconstructed using the GAL 
or the SRTD. When both GAL and SRTD information are available, the SRTD reconstruction 
is used. 

• GAL - The position of a shower in the GAL can be determined using the pulse-height 
sharing between left and right PMTs and between neighbouring cells. The position reso- 
lution on the face of the GAL is about 1 cm. The position measurement was checked by 
matching tracks and positron candidates. The resolution of the track-GAL matching is 
found to be 1.1 cm and systematic biases are less than 2 mm (see Fig. |I]a and b). 

• SRTD - The position in the SRTD is determined from the centres of gravity of the pulse- 
heights of the strips in the X and Y planes of the SRTD. The measured SRTD position 
resolution is 0.3 cm. For a subset of these positrons the position measured by the SRTD 
can be compared to the result from the track measurement. The resolution of the track- 
SRTD matching is measured to be 0.5 cm and systematic deviations are less than 0.1 cm 
(see Fig. |l|c and d). 

The coordinates of the event vertex are determined from tracks reconstructed with the GTD 
and VXD. The Z coordinate {Z^ertex) is determined on an event-by-event basis. Since the 
transverse sizes of the beams are smaller than the resolutions for the X and Y coordinates of 
the vertex, the beam X and Y positions averaged over a single beam-fill are used. For events 
which do not have a tracking vertex, the Z coordinate is set to the nominal position of the 
interaction point {Z = +6 cm). 

The probability of finding a vertex depends on the number of particles which traverse the 
tracking detectors and thus mainly on the hadronic angle 7^. Figure |I|e, which shows the 
fraction of events having a reconstructed vertex as a function of 7^, exhibits the expected 
behaviour. For 7^ > 70° the vertex finding efficiency is greater than 95%. It decreases to 
around 50% for events with 7^ ~ 20°. The MG simulation reproduces the observed behaviour 
reasonably well. Relative differences of the order of 5-10% are observed in the region below 
7^ = 30°. The drop in vertex finding efficiency in the region around 7^ = 174° results from 
the combination of the characteristics of the DIS final states and the acceptance of the GTD 
in the rear direction. The effect is not fully reproduced by the MG simulation. All differences 

^ Contrary to the previous analysis |^] , in this analysis we allow the cut on the probability to vary with the 
calorimeter energy associated with the positron. 
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between data and MC simulation are taken into account in the determination of the systematic 
uncertainties on the final results. 

Figure ^ shows the Z resolution of the vertex reconstruction as a function of 7^. For 7^ > 40° 
the resolution is a^^"''^^ = 0.2 cm, while for events with 7^ close to zero f\ it becomes 
track _ 2.5 ± 0.4 cm. For events where no vertex is reconstructed, the effective resolu- 
tion is determined by the length of the interaction region, which is a = 12 cm. It is also 
possible to determine the Z coordinate of the vertex from the measurements of the arrival times 
of energy deposits in the FCAL In this measurement the resolution is = 17 cm, and 

it is only used for studies of systematic errors. 

For polar angles 45° < 9^ < 135° and positron transverse momenta pTe > 5 GeV the tracking 
efficiency is greater than 98%. Thus when the positron angle, as determined by the GAL, is 
smaller than 135° a track is required to match the positron identified in the CAL and the polar 
angle is taken from the track parameters. A successful match is obtained when the distance 
between the extrapolated impact point of the track on the face of the CAL and the position 
determined by the CAL is less than 5 cm. In this region the matching efficiency is above 97% 
and the combined efficiency for tracking and matching is greater than 95%. 

The resolution of the scattering angle measurement is cf/*™ = 2.0 mrad for positrons recon- 
structed in the SRTD, cr£"^^ ^ 6.8 mrad for positrons determined by the CAL and Ug™'^'^ ~ 
3.4 mrad for positrons with a matched track. For events without a tracking vertex there is an 
additional error of ci^^^*^^ ^ 80 (sin 6e) mrad. 



4.3 Positron energy determination 

The scattered positron loses energy in the passive material in front of the CAL. In the region 
relevant for this analysis this passive material constitutes about 1.5 radiation lengths except in 
areas around the rear beampipe, 9 > 170°, and the solenoid support structure, 130° < 9 < 145°, 
where it is up to 2.5 radiation lengths. In the analysis, the measurement of the scattered 
positron energy by the CAL is corrected for the energy loss in the passive material. 

The correction for the CAL measurement of the scattered positrons can be determined directly 
from the data using the following subsamples. 

• At low y, the scattered positron energy is kinematically constrained to be close to the 
positron beam energy and to be primarily a function of the scattering angle. For these 
events, called Kinematic Peak (KP) events, the mean positron energy is determined from 
the scattering angle {E'^ ^ {i-mle ) ) within 0.5%. The KP events are selected by 
requiring = ^ < 0.04, and provide a calibration at E'^ ^ Ee = 27.5 GeV for 
9 > 135°. 

• For QED Compton events {ep — >■ e'yp) observed in the main detector, energies of the 
positron and the photon can be predicted precisely from the measurement of their scat- 
tering angles since the transverse momentum of the scattered proton is small. QED 
Compton events provide a calibration at 5 < -E^ < 20 GeV and 9 > 160°. 

■^The lines of constant 7^ in the (a:, Q^)-plane are given in Fig. |^ for 7^ = 40°, 90° and 135°. 
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• In events from DIS p° production, ep epp° (p° tt+tt^), the angle of the positron and 
the momenta of the 7r+ and vr", as measured with the CTD, give a precise determination 
of the positron energy. DIS p° events provide a calibration at 20 < E'^ < 25 GeV and 
e > 160°. 

• In the range 30° < 6 < 150°, momenta of the positrons can be determined by the CTD. 
Although the CAL energy resolution at E'^ > 10 GeV is superior to the momentum reso- 
lution of the CTD, the track momenta averaged over several events give an independent 
check on the energy measurement of the CAL. 

In the fiducial volume of the SRTD, 9 > 167°, the correlation between the energy lost in the 
passive material in front of the calorimeter and the energy deposited in the SRTD is used to 
correct the calorimeter energy measurement (SRTD method). The corrections are determined 
using the QED Compton and the KP samples. The DIS p° events are used as checks on the 
correction. This procedure is described in detail elsewhere 

In the RCAL outside the SRTD fiducial volume, 130° < 9 < 167°, the observed energy shift of 
KP positrons from the value expected from kinematics, combined with test beam results of the 
energy loss of electrons in passive material, are used to determine the correction (KP method). 
This procedure is described in detail elsewhere [Q. Using the data from regions where both 
the KP method and the SRTD method can be applied, the uncertainty in the KP method of 
correction is found to be the same as that of the SRTD method. 

The uncertainty in the energy determination in the RCAL after these corrections is 2% at 10 
GeV linearly decreasing to 1% at 27.5 GeV. 

In the region covered by BCAL and FCAL, 9 < 130°, comparison with the measured mean 
CTD momentum is used to correct the CAL energy measurement (track method). The ener- 
gies corrected by the track method have an uncertainty of 2% independent of the energy of 
the positron. Where the KP method and the track method corrections can both be applied, 
consistent results are found. 

The resolution of the positron energy measurement can be determined from 
the data using QED Compton, DIS p°, and KP events, and is found to be in the range 
a/E = (20 — 27)%/-y/i? (GeV), depending on the trajectory of the positron through the passive 
material in front of the CAL. 

Studies identical to those described above for the data were also performed with MC event 
samples. The resulting corrections were applied to the DIS MC sample. Thus, any difference 
in the CAL response between MC simulation and data is absorbed in the corrections. The 
resolution of the positron energy measurements in the MC simulation is adjusted to reproduce 
the results found from the data. 

The separate positron energy corrections for MC simulation and data allow the determination of 
relative calibrations for the CAL in MC simulation and data. Before determining the hadronic 
energy, adjustments, based on the difference in the magnitude, position dependence and energy 
dependence of MC and data corrections, are made to bring the energy scale for data and MC 
simulation into agreement. 
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4.4 Hadronic energy determination 



The hadronic energies, unhke the case of the positron energy, are not corrected for energy loss in 
the passive materiaL Instead, the transverse momentum of the positron, p^ei calculated using 
the positron energy corrected as described in the last section, is compared to the pxh of the 
hadrons in both the detector simulation and data. From this comparison, uncertainties in the 
determination of the hadronic energy are estimated. The mean as a function of 7^ agrees 
within 3% between MC simulation and data for the entire range of kinematics covered in this 
paper. The distributions in bins of 7^:^ are compared (see Sect. ^ for MC simulation and 
data, and are in good agreement. An uncertainty of ±3% is assigned to the hadronic energy 
measurement, based on these comparisons. 

5 Kinematic reconstruction 

In deep inelastic scattering, e{k) + p{P) + X, the proton structure functions are 

expressed in terms of the negative of the four-momentum transfer squared, Q^, and Bjorken x. 
In the absence of QED radiation. 



where k and P are the four-momenta of the incoming particles and k' is the four-momentum 
of the scattered lepton. The fractional energy transferred to the proton in its rest frame is 
y = Q'^/{sx) where s is the square of the total center of mass energy of the lepton-proton 
collision {s = 90200 GeV^). 

The ZEUS detector measures both the scattered positron and the hadronic system. The four 
independent measured quantities E'^, 9e, Sh and pxh, as described in the previous section, over- 
constrain the kinematic variables x and (or equivalently, y and Q"^). 

In order to optimise the reconstruction of the kinematic variables, both the resolution and 
robustness against possible systematic shifts (stability) of each measured quantity must be 
considered. 

For the present analysis, a new method (PT) is used to reconstruct the kinematic variables. 
The PT method achieves both superior resolution and stability in x and in the full kine- 
matic range covered, in comparison with reconstruction methods used in our previous structure 
function measurements. 

5.1 Characteristics of standard reconstruction methods 

As discussed in Sect. ^ the positron variables, E'^ and 9e are measured with high precision, and 
the systematic uncertainties are small. The kinematic variables calculated from these quantities 
are given by: 



Q' = -q' = -{k 
X = 




(6) 



(5) 



Ve = I- 




(7) 
(8) 
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This method of reconstruction, the electron (EL) method, gives good results at high where 
E'^ is significantly different from the positron beam energy E^,, but at low y {y ^ 0.1) both the 
resolution and stability of ye become poor. 



The Jacquet-Blondel ( JB) method of kinematic reconstruction only uses information from 



the hadronic energy flow of the event: 



2Ee 



(9) 



Qls = (10) 

Qe has better resolution than over the entire kinematic range while at low y {y ^ 0.04), 
yjg has superior resolution in comparison with ye. However, unlike the case of the positron 
measurement, the energy loss of the hadronic system in the passive material in front of the 
CAL and the loss through the beam holes have to be determined using a simulation of the DIS 
flnal states and of the detector effects. This introduces substantial systematic uncertainties. 



The Double-Angle (DA) method ||28| combines the information from the scattered lepton with 



that from the produced hadronic system: 

n2 = ,p2 sin7^(l + cosge) ... ^ 

' sin 7^ + sin Oe - sin(7^ + Oe) ' ^ ^ 

singe(l-cos7„) 
sin 7^ + sin Oe - sin(7^ + ^e) ' 

In this method the hadronic measurement enters the kinematic reconstruction through the 
variable 7^ (see eq. §), which depends on the ratio of the measured quantities 5h and pxh- As 
a result, uncertainties in the hadronic energy measurement tend to cancel leading to a good 
stability of the reconstructed kinematic variables. At the highest ?/, however, ye has better 
resolution than y^^. At low y ( < 0.04), and low Q^, the noise in the CAL, which affects both 
PTh and 5h has a large effect on and j/^^. This region was avoided in our previous analyses, 
which used the DA method, by imposing the cut y > 0.04. 



5.2 PT method 

The PT method provides an improved measurement of y and by an efficacious combination 
of the information from the measurements of both the hadronic system and the positron. 

This is done in two steps. In the first step, the transverse momentum balance between the 
positron and the hadron system is used event by event to correct the hadronic energy measure- 
ment of y by using a functional form derived from the MC simulation. In the second step, the 
advantages of the EL, JB and DA methods are combined into a single reconstruction method 
for the full kinematic range. 
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5.2.1 The measurement of y at low y 



In an ideal detector, transverse momentum conservation requires pxh = Ptc up to the negligible 
transverse momentum carried by QED initial state radiation. Since the positron is measured in 
the calorimeter, energy carried by QED final state radiation is mostly included in the measured 
positron energy. Since pre is a precisely measured quantity (see above), prh — Pre gives, for 
each event, the difference between the measured and the true hadronic transverse momentum. 
The difference can result from hadronic energy loss through the beam hole and/or in the 
passive material in front of CAL, as well as from the finite energy resolution of CAL. In 
first approximation the hadronic energy mismeasurement affects both pth and 6h in the same 
way. Indeed, if the hadronic system consisted of the proton remnant carrying no transverse 
momentum and a well coUimated current jet, then an accurate measurement of y would be 
obtained by: 

yii)=yjB/i'^)- (13) 

In this idealised case, a precise measurement of y is achieved by determining the hadronic 
energy correction directly from the data (on an event by event basis) rather than by a MC 
simulation. 

In standard DIS events, in addition to a current jet, there is hadronic energy fiow in the region 



between the current jet and the direction of the proton remnant [1/]. As a result, the simple 



relationship given by eq. |T3| no longer holds and is replaced by: 

= Vjs/C (14) 

where C is a correction function. In the PT method, C is determined from MC simulation of 
final states convoluted with the detector simulation, as a function of the measured quantities 

1^, pTh and 7^ giving C{^,pTh,lH)- 

The ratio averaged over pxh, {^^), as obtained from the MC simulation, is shown in Fig. El 

Vgen ^ Vgen, ' — ' 

as a function of ygen is the generated value of y. For small values of 7^, (^^) rises almost 
linearly with This behaviour is as expected from eq. p!3| and follows from the fact that the 
hadronic energy deposit comes mainly from the current jet when 7^ is small. Note that the 
quantities pTh and yjg are rather insensitive to the energy loss in the forward beam hole. The 
observed positive offset from unity of (^^) at — = 1 results from the CAL noise. R 

As 7^ increases (^^) falls below unity mainly due to energy loss in the passive material in the 
detector. With increasing 7^ an increasing fraction of prh is carried by particles produced in 
the region between the current jet and the proton remnant. The energy loss of these particles 
enters directly into the measurement of pTh but to a lesser extent into the measurement of yjg 
(see eqs. |^ and ||). This leads to a less linear dependence of (^^) on 

As PTh increases the energy fiow between the current jet and the remnant becomes relatively 
less in comparison to the energy fiow in the current jet. Thus, as a function of increasing pth, 
the dependence of becomes more linear (not shown in figure). 

PTe ^ Vgen ^ O / 

* As an example of the effect of calorimeter noise on ( ) at low , consider the measurement of an 

event with true y = 0.01 and true — 10 GeV^. If one calorimeter cell in the RCAL produces a noise 
signal of 100 MeV, the contribution to the measured tjjg amounts roughly to '^^^q^^ = 0.004. Ignoring the 
effects of energy loss and smearing, is reconstructed 40% larger than the true y. The contribution of the 
same noise signal to the measured pxh, on the other hand, is small (< 0.1 GeV) in comparison with the true 
PT = Qy(r=y) = 3.15 GeV. 
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The correction function C{^,pTh,lH) parameterised using second order polynomials to 
describe the MC simulation of in the range 0.5 < — < 1.5 which contains approximately 
90% of the events. The curves in Fig. ^ show C, averaged over p^h, as a function of for 
several values of ■y^ . 



Application of the correction function, eq. |Tj, determines y based on the measitred values of 



PTh and 7^, which reduces the dependence of the y measurement on the MC simulation. At 
the same time, by compensating for the deviations of the hadronic energy measurement from 
the true values for each event individually, the resolution of the y measurement is improved. It 
will be shown in Sect. [7?T| , in the studies of systematic uncertainties, that the structure function 
results do not depend strongly on the details of the MC simulation of the DIS final states used 
in determining the correction function. 



5.2.2 Kinematic measurement in the full y range 

When the current jet points in the backwards direction (7^ > 90°), corresponding to the region 
of high y, the hadronic system becomes less collimated, and the particle loss through the rear 
beam hole is not negligible. These particles carry a small amount of transverse momentum but 
make a large contribution to the decrease in 6h and therefore to the decrease in yjg. As a result 
the ratio is less effective for determining the correction to the hadronic measurement of v. 
This can be seen in Fig. ^ where for 7^ > 120°, there is little dependence of (^^) on 

The measurement of y from the positron (see eq. |^ at high y does not suffer from these 
deficiencies and provides an accurate measurement. Applying the S-correction, first introduced 



in ref. p9[ and used by HI [g], and ye are combined: 



1/(2) = 1/(1) — n • (^^) 

At this stage, the best estimate for the transverse momentum of the hadronic system is pTe 
and for the hadronic sum J2h{Eh — Pzh) is 2E'^y{2)- This allows the calculation of the hadronic 
angle, 7^,,, as follows: 

vie - ^Elyl) 

^QS7PT= 2 , 4Z.2^2 • (16) 
Pre + ^-^el/(2) 

Since the Double-Angle method is the method least sensitive to energy scales, the kinematic 



variables are calculated using eqs. 11 and 12 substituting 7^^, for 7 



The kinematic variables determined in this way will be denoted by the subscript PT. Figure ^ 
shows the distributions of and for several ranges in 7 . By construction, the 

Vgen ygen Vgen Vgen 

correction from y to centers the distribution around unity, while the further corrections 
yield improvements in resolution. Figure |3|b shows the distributions of together with 

' — ' Vgen Vgen 

and The corresponding distributions for the reconstruction are shown in Fig. |3|c. Both 
in y and the resolution of the PT method is superior to all other methods. Also shown in 

2 

Fig. ^ is the result for Q(2)(~ iz^^)- The resolution of this alternative method to reconstruct 
is inferior to the resolution of the PT method. 
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For large y one can check the method by comparing with Uf.. Figure ^ and b show the 
distributions of for the regions 140° < 7^ < 160° and 160° < 7^ < 180° where Ue is 
expected to give a reasonable estimate of y . Also shown in Fig. ^ are the expectations from 
the MC simulation. The r.m.s. widths (a = 9% for 160° < 7^ < 180° and a = 20% for 
140° < 7^ < 160°) of the distributions are in good agreement with the MC simulation and 
are dominated by the error on the measurement of ye ■ There are small shifts of about 2% 
between the data and MC simulation. These shifts are included in the study of the systematic 
uncertainties. 



6 Event selection 



The following cuts were used to select NC DIS events: 



A positron candidate identified as described in Sect. |41 . 



-Eg > 10 GeV, where E'^ is the corrected positron energy. This cut ensures high and well 
understood positron finding efficiency and suppresses background from photoproduction. 

38 GeV < 5 < 65 GeV, where 6 = J2i{Ei — pzi)- Here the sum runs over all calorimeter 
cells, i.e. including those belonging to the identified positron. This cut removes events 
with large initial state radiation and further reduces the background from photoproduc- 
tion. 

A track match for 6 < 135°. This condition suppresses events from cosmic rays, halo- 
muons, photoproduction and DIS events where an electromagnetic shower has been falsely 
identified as the scattered positron. 

ye < 0.95. This condition removes events where fake positrons are found in the FCAL. 



A total of 680283 events pass the above cuts. 



• Box cut. Events with a scattered positron impact point in the RCAL inside a box of 
26 cm X 26 cm around the beampipe are rejected. This ensures that the impact point is 
at least 2.5 cm away from the edge of the RCAL and therefore guarantees full containment 
of the electromagnetic shower in the calorimeter. 

• —50 cm < Zyertex < 100 cm. This cut is performed for events with a reconstructed 
tracking vertex and suppresses beam-gas background events and the small fraction of 
the events where the vertex position is incorrectly measured. Events without a tracking 
vertex are assigned the mean vertex Z position and are accepted. 



After these additional cuts a sample of 443421 events remains. 



• 1^ > 0.3 — 0.0017^ (7^ in degrees). For the PT method to yield a reliable measurement 
of X and the loss in hadronic transverse momentum must be limited. There are events 
where a substantial part of the current jet remains in the forward beampipe. These events, 
which are produced at small 7^ , can be falsely reconstructed at large 7^ due to the CAL 
noise. This cut suppresses such events. 
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A total of 400627 events pass all the above selection cuts. For the accepted events Fig. ^,b 
show, respectively, the measured energy and scattering angle of the positron. Figure ^ shows 
the reconstructed vertex distribution before applying the vertex cut and Fig. |^-f show the 
distribution of ^ , where the ^ cut was removed, for three ranges in 7^ . The energy, angle, 
vertex and transverse momentum distributions are well reproduced by the MC simulation (also 
indicated in the figures). To obtain these distributions the MC events have been reweighted to 
the structure function obtained from the QCD NLO fit described below. 

Figure 10 shows the distributions and 7^^, together with the MC distributions 

normalised to the integrated luminosity of the data. Good agreement is obtained except for 
the region around y ~ 10~^ where the MC undershoots the data by up to ~ 5%. Equivalent 
deviations are observed in the Xp^. and 7pj, distributions. This difference is not concentrated 
at a specific value of and has a negligible effect on the values of the extracted structure 
function. 

Figure |^ shows the distribution of events in the plane. The {x^Q"^) bins used for the 

determination of the structure function are shown in Fig. 0b- The bins have been chosen 
commensurate with the resolutions. At large and also at low y larger bin sizes have been 
chosen to obtain adequate statistics in each bin and to minimise bin-to-bin migrations as a result 
of the non-Gaussian tails. Furthermore, large acceptance {A > 20%) and purity (P > 30%) Q 
are required for each bin. The lowest values of acceptance occur at the lowest values where 
the box cut becomes effective. The bins with lowest purity occur at the lowest y values. In the 
majority of bins the acceptance is greater than 80% and the purity is greater than 50%. 

Figure § shows the distributions in ( ^PT-xgen \ ^QptJ^msil\ for a number of bins in the (x,Q^) 
plane, where Xgen and Q^en ^^re the generated values of x and Q^. The resolutions vary smoothly 
with X and Q^. The resolution in x decreases from ~ 25% for low y {y = — ) to cr^. ~ 10% 
for large y. The resolution in is aQ2 ^10% for large y and improves to aQ2 ^5% for 
intermediate and low y. At the lowest y values tails appear due to the vertex finding efficiency 
and resolution. 

The numbers of events in each of the bins are given in Table |I|. 

The final sample contains a small number of background events which are not due to deep 
inelastic neutral current scattering: 

• Non-ep background. The level of background not associated with ep collisions is deter- 
mined from the number of events observed in unpaired or empty bunches. This back- 
ground is subtracted statistically taking into account the appropriate ratios of bunch 
currents and numbers of bunches. It amounts to less than 1% in all bins. For events in 
the bins of > 1000 GeV^ all events were scanned visually and no non-ep events were 
found. 

• Photoproduction background. The events from the photoproduction MC code PYTHIA 
were analysed in the same way as the data and the number of events passing the selection 
cuts in each bin was determined. This number is subtracted taking into account the ratio 

^Acceptance is defined as the ratio of the number of events which are reconstructed to the number of 
generated events in a bin. Purity is defined as the fraction of reconstructed events in a bin which originated 
from the same bin. 
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of equivalent luminosity of the data and the MC samples. As a check the S distributions 
for each of the bins were fitted to a shape extracted from MC DIS events and a Gaussian 
shape, which parameterises the photoproduction events (the procedure is described in [|1|). 
The 6 distribution for y > 0.2 is shown in Fig. Also given are the MC expectations 
for DIS (solid histogram) and for photoproduction (dashed histogram). The flattening 
of the MC prediction for photoproduction at 5 < 38 GeV is artificial and results from 
the cut of W > 190 GeV applied at the generator level. In the region 5 < 38 GeV 
the photoproduction contribution can be measured directly through events tagged by the 
LUMI positron tagger (the acceptance for these events tends to zero for 5 > 40 GeV). 
The result of this measurement is shown in Fig. ^. Also shown is the result of the fit 
to the data (dashed-dotted curve) using the shape from the DIS MC distribution plus a 
Gaussian contribution (dotted curve) for the photoproduction background. The Gaussian 
contribution agrees with the MC photoproduction estimate for (5 > 38 GeV and with the 
measured photoproduction background. For > 1000 GeV^, the visual scan found at 
most two events showing topologies consistent with photoproduction events. 

The total background contribution to each bin is given in Table ^ The maximum background 
fraction of about 6% occurs in bins of high y. 



7 Proton structure function F2 

In deep inelastic scattering the double differential cross section for inclusive e~^p scattering is 
given in terms of the structure functions Ff. 

\Y^F2{x,Q')-y'F^{x,Q')-Y.xF,ix,Q')] (1 + 5,(x, Q^)) (17) 



dx dQ^ xQ^ 

where Y± = 1 ± {1 — y)^ and x and are defined at the hadronic vertex . In this equation 
F^ is the longitudinal structure function, F3 is the parity violating term arising from the Z° 
exchange and 6r is the electroweak radiative correction. Since all of these contributions to the 
cross section are expected to be small in the kinematic region of the present measurement, this 
can be rewritten as: 

The F2 structure function itself contains contributions from virtual photon and Z° exchange: 

where Mz is the mass of the Z° and F2'^, F^*^ and ^2"* are the contributions to F2 due to 
photon exchange, Z° exchange and 'yZ° interference respectively. We finally write: 

^^"'^"-Fr(l + 5.)(1 -S,- 6s){l + 6r). (20) 



dx dQ' xQ^ 
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The corrections, ^r,z,3,i,, are functions of x and but are, to a good approximation, indepen- 
dent of F2, i.e. they are insensitive to the parton density distributions. 

In this analysis we determine the structure function F2. In the Standard Model the difference 
between F2 and Fl^"* is expected to be less than 1% for values of < 1000 GeV^. At higher 
the difference becomes progressively larger. In this region therefore we present the values of 
¥2'^ as well as F2. These are extracted from F2 through the correction 5z using the MRSA[p[ 
parton distribution functions. 

An iterative procedure is used to extract the structure function F2. Monte Carlo events are 



generated according to eq. and so contain transverse and longitudinal photon and Z° contri- 
butions as well as radiative effects. In each step of the iteration the value of F2 in a bin can be 
obtained from the ratio of observed events in the data to the number of events observed in the 
MC simulation and the calculated values of 5r,z,^,L- ■'■^ ^^^^ ^^^P acceptance correction is 
taken from the MC simulation which uses the parton distributions given by MRSA. The result 
for F2 from this first iteration is used for a QCD fit using the DGLAP |Q evolution equations 



in next-to-leading order (QCD NLO fit) in a manner very similar to that described in ^T[|. The 
evolution uses massless quarks of three flavours in the proton and the charm quark coefficient 
functions from references WA, RH] to ensure a smooth crossing of the charm threshold. The 



NMC pi data for > 4 GeV^ are used to constrain the QCD NLO fit at high x ( > lO'^^ 



The fit also includes the low Q 1994 data from this experiment P]. 

The parton distributions from the QCD NLO fit are used to recalculate the cross section, 
including the new Fi contribution as a function of x and Q^. The MC events are reweighted 
using the calculated cross sections and the true hadronic x and Q^. This reweighted MC is then 
used to perform the acceptance correction bin by bin to the data, leading to a new estimate of 
F2, specified at the mean values of x and in each bin. The procedure is repeated until the F2 
values from two consecutive iterations change by less than 0.5%. The final result is reached in 
three iterations. The functional form of the QCD NLO fit is used to quote the results for F2 at x 
and values which are chosen to be convenient for comparisons with other F2 measurements. 
Using alternative parametrisations for this correction has negligible effect on the measured F2. 
It should be noted that the QCD NLO fit is used here only as a parameterisation to obtain a 
stable acceptance correction and not to perform a detailed QCD analysis. 

The statistical errors of the F2 values are calculated from the number of events measured 
in the bins, including those from the background subtraction, and the statistical error from 
the MC simulation. Since we are using bin by bin corrections, the correlations of statistical 
errors between the F2 measurements enter only via the finite statistics of the MC sample. The 
correlations are small given the relatively large MC sample used in this analysis. A correlation 
between the F2 values of neighbouring bins is present due to acceptance and smearing effects. 
The sensitivity of the measured F2 to these effects has been checked by comparing the F2 
obtained from the first iteration with the final value. The changes were within the statistical 
errors of the final F2 values. 



The correction 5^ is determined using the QCD prescription for F^ ||35| and the parton distri- 
butions from the QCD NLO fit. It is significant only in the region of large y where it reaches 
12%(see Table p. The 5^ and ^3 corrections are calculated from the MRSA parametrisation 
of the parton distributions. These corrections are negligible for values of below 1000 GeV^. 
For values of above 1000 GeV^ they strongly increase with increasing and have a slight 
dependence on y. For the two highest bins the corrections are: 
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g2 = 3000 GeV^, 

X = 0.20: 
X = 0.08: 6^ 

2 - 5000 GeV^ 

X = 0.20: 5^ 
X = 0.08: 6^ 



Q 



+4% and ^3 
+5% and ^3 



+5% 
+8% 



-7.5% and ^3 = +11% 
-8% and ^3 = +22%. 



The F2 values are corrected for higher order QED radiative effects not included in HERACLES. 
These corrections, including soft photon exponentiation, were evaluated using the program 
HECTOR in the leading log approximation They vary smoothly with between 0.2% 
and 0.5%. 



7.1 Systematic uncertainties 

Several factors contribute to the systematic uncertainties in the measurements of F2. In the 
following they are subdivided into six independent categories. For each category we summarise 



the checks that have been performed to estimate the size of the uncertainty. Figure |T0| shows 
the size of the systematic uncertainties plotted as a function of y for the six categories and also 
the total systematic error. For reference (see appendix) each test apphcable at < 100 GeV^ 
OT y < 0.01 is numbered; the number is given in brackets {}. 



1 Positron finding and efficiency 

• The positron identification efficiency is varied within the uncertainty found in the 
QED Compton study which is ~3% for positron energies E'^ = 10 GeV and negligible 
for energies above 14 GeV {1,2}. The effects on F2 are negligible except in the lowest 
X bins where they reach ~3%. 

The parameters of SINISTRA94 are varied to check the effects from overlaps of 
energy deposits from the hadronic final state with those of the scattered positron 
{1}. This results in changes of F2 of ~4%, again in the lowest x bins. 
For Of. < 135°, where a track is required to match with the CAL positron position, 
a comparison of the number of rejected events in MC and data can accommodate a 
possible mismatch in efficiency of at most ±2.5%. 

In addition, an alternative electron finder employed in a previous analysis is used 
Consistent results are obtained over the full energy range. 

2 Positron scattering angle 

• Changing the box cut from 26 cm x 26 cm to 28 cm x 28 cm in both data and MC 
{3,4} has a small effect on the F2 at low values of Q'^{< 15 GeV) and are negligible 
elsewhere. 

The systematic uncertainty in the position of the scattered positron is estimated 
from the difference between the extrapolated track position at the face of the CAL 
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and the position found from the CAL or the SRTD. In the region of the SRTD 
these uncertainties are about 1 mm. They are 2 mm outside the SRTD. The largest 
systematic error is obtained if the relative distance between the two halves of the 
RCAL is changed by ±2 mm {5,6}. 

For positions determined by the RCAL we additionally changed the distance of the 
measured impact point with respect to the beam-line by ±2 mm {7,8}. 
For positron angles measured by the tracking detector, the scattering angle Of, is 
changed by the estimated systematic uncertainty of the track angle of ±0.2°. 
These changes were made in the data while leaving the MC simulation unchanged. 
This had an effect of about 1% on the measured F2. 

• The absolute vertex reconstruction efficiency in the MC simulation is decreased by 
1% overall {9} or increased by 3% for 7^ < 40° {10}. These changes are motivated 
by the study of the efficiency as a function of the hadronic angle (see Fig. |l]e). Also, 
the cut on the Zyertex position is tightened by demanding -28 cm < Zyertex < 40cm 
to estimate a possible uncertainty from the luminosity of the events in the satellite 
bunch {11}. Effects on the measured F2 are typically 1-2%. 

As an additional check, we required that the vertex determined from the CAL timing 
was -40 cm < ^^erteaf < 30 cm for the events that had no tracking vertex. The effects 
on the measured F2 are similar to those described above and are not included in the 
final systematic uncertainties. 

3 Positron energy scale 

• The systematic uncertainties in F2 due to the uncertainties in the absolute calorime- 
ter energy calibration for scattered positrons is estimated by changing the energy 
scale in the RCAL in the MC simulation by 2% at 10 GeV linearly decreasing to 1% 
at 27.5 GeV. For BCAL the energy scale is modified by ±2%. The magnitude of these 
shifts represents our present understanding of the positron energy measurement (see 
Sect. ^) {12,13}. Systematic variations in F2 of 1-2% are observed. 

4 Hadronic energy measurement 

• The energy scale in the MC simulation is changed by ±3% while leaving the positron 
energy scale unchanged {14,15}. The value of 3% is based on detailed comparisons 
of the distributions of the quantity from data and MC (see Fig. |]d-f ) and on the 
deviation of the mean of from unity for high y (see Fig. The effects on F2 
are typically 4% in the bins at very high and very low while the effects are about 
1% at moderate y. In addition the energy scales of the different calorimeter parts 
are changed independently by 3% in turn {16-19}, which resulted in similar effects 
on F2. 

• The noise cut on isolated EMC(HAC) cells of 100(150) MeV was changed in the MC 
simulation by ±10 MeV which reflects the uncertainty in the MC simulation of low 
energy deposits in the CAL and the modeling of the uranium radioactivity {20,21}. 
Variations in F2 of up to 25% are observed in the bins at very low y. Elsewhere the 
variations in F2 are within 2%. 
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5 Hadronic energy flow 



The correction function, C{^^,pTh,lH)^ depends on tlie liadronic energy flow. The 
correction function was also determined using a MC simulation which uses the 
LEPTO MEPS model for the description of the hadronic final state, instead of 
CDMBGF. This model also provides a reasonable description of the hadronic energy 
flow in DIS events [0. The new correction function was applied to the data. Ac- 
ceptance corrections were then made with the sample of CDMBGF MC events with 
the original correction function applied {22}. The changes in F2 are typically 2% in 
all bins. 

The hadronic energy flow for diffractive events is different to that for non-diffractive 
events. To investigate the sensitivity of the PT method to the size of the diffrac- 
tive component we reweighted the diffractive scattering cross section in the MC 



simulation such that it agrees with the measurements of |^ while leaving the cor- 
rection function C unchanged {23}. As a separate check, to allow for possible addi- 
tional diffractive- like events we increased the diffractive cross section by an additional 
100%{24}. Effects at the 2% level in F2 are observed in each check. 

• The fraction of events removed by the cut is sensitive to the amount of pTh lost 
in the forward beampipe and so also to the details of the fragmentation. We removed 
the 1^ cut, while extending the correction curves to include the region of very low 
|^{25}. This resulted in changes of 3% percent at low y. 

• We determined the correction functions without allowing for a dependence on pTh 
{26}, resulting in changes of less than 2% in F2 . Reducing, by a factor of two, 
the size of 7^ intervals in which the correction function is determined resulted in 
changes of typically l-2%{27}. 

6 Background subtraction 

• As described above, the photoproduction background as given by the PYTHIA MC 
simulation agrees with the determination from the fitting procedure described in 
section ^ (see Fig. H), which in turn agrees with the photoproduction contribution 
extracted from events tagged by the LUMI positron detector. Uncertainties in the 
LUMI positron detection efficiency, the assumed shape in the fitting procedure and 
the positron misidentification in the PYTHIA MC simulation lead to an overall 
uncertainty in the photoproduction background of ±50%{28,29}. This results in an 
uncertainty in F2 which is at most 3% for high y and negligible elsewhere. 

The total systematic uncertainty for all bins is determined by adding in quadrature separately 
the positive and negative deviations from the above mentioned categories. 

At high (> 100 GeV^), a significant contribution to the apparent systematic deviations 
is due to the limited number of events in the data. From the lower region we know that 
the systematic uncertainties depend mainly on y. Therefore, all bins with > 100 GeV^ 
and y > 0.01 are combined, for the systematic error determination, in fixed y intervals. The 
uncertainty in the track matching efficiency contributes only at > 350 GeV^ and is taken 
into account separately. 

Figure ffUl shows the positive and negative systematic errors for each bin with < 100 GeV^ as 



a function of y. For the majority of the bins the total systematic uncertainty is below 5%. As 



18 



y increases above y = 0.5 the uncertainty grows to around 10%. For y < 0.01 the uncertainty 
increases mainly because of the CAL noise. 

The combined results for the bins with > 100 GeV^ are also shown. 

The systematic errors do not include the uncertainty in the measurement of the integrated 
luminosity (±1.5%), the overall trigger efficiency (±1.0%) or the uncertainty due to higher 
order electroweak radiative corrections (±0.5%). These effects lead to a combined uncertainty 
of 2% on the overall normalisation of F2. 

In addition to the above studies a complete analysis was done using the electron (EL) method 
of kinematic reconstruction. In the region of high y, where this method is reliable, the results 
were compatible with this analysis. 



8 Results 



The values of F2 are given in Tables and ^ together with their systematic and statistical 
uncertainties. In the appendix we give the F2 values with < 100 GeV^ or y < 0.01 obtained 
from each systematic check used in the calculation of the systematic uncertainties. Also shown 
in Tables |I] and ^ are the number of events and the estimated number of background events in 
each bin as well as the correction due to the longitudinal structure function, 5^ . The F2 values 
are displayed versus x for fixed values of in Fig. 11 together with the low data from our 



experiment and the high x results from NMC [0] . The rise of F2 for x — > is measured with 
much improved precision. Also shown is the result of the QCD NLO fit used for the acceptance 
determination. The data from this analysis at large x reach the x range covered by the fixed 
target experiments. In the overlap region the agreement is good. The results agree well with 
the recently published results of the HI collaboration |^ in the kinematic region where the data 
sets overlap. 

Figure shows the F2 values as a function of for fixed x. Scaling violations are observed, 
which decrease as x increases. The Q^-range of the ZEUS data has increased substantially with 
respect to our previous data. The data now span more than a decade in Q"^ for x values as low 
as X = 4 • 10~^ and almost three decades at x ~ 0.1. 



In Fig. 13 the data are shown together with the NLO predictions of GRV94 |38| and the 



NLO parameterisations MRSA' [^] and CTEQ3 The data of the fixed target experi 



ments [0, ^ ^ are also shown. At large all parameterisations represent the data well. 



For < 45 GeV^ CTEQ3 undershoots the measured F2 values in the region 10 ^<x<10 ^. 
MRSA' overshoots the data below x = 10"^ at low (Q^ < 8.5 GeV^), whereas GRV94 over- 
shoots the data at < 18 GeV^ and x < 10~^. Our QCD NLO fit reproduces the data over the 
full kinematic range indicating that NLO DGLAP evolution can give a consistent description 
of the data. The details of this fit will be included in a forthcoming paper. 



9 Conclusions 



We have presented new measurements of the proton structure function F2 from an analysis of 
inelastic positron proton neutral current scattering using data obtained with the ZEUS detector 
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at HERA during 1994. A new method for determining the event kinematics has allowed us to 

measure F2 over a substantially larger phase space than before by this experiment and with 
systematic uncertainties reduced to below 5% in most of the (x, Q^) space. The data cover Q"^ 
values between 3.2 and 15000 GeV^ and x values between 5 • 10~^ and 0.8. 

At large x values and Q"^ values up to 70 GeV^, where the new data reach the x-range covered 
by fixed-target experiments, good agreement with these experiments is observed. The data 
show the rise of F2 towards small x with much improved precision. Strong scaling violations 
are observed for x < 0.02. The measured x-Q^ behaviour of F2 can be described by QCD using 
NLO DGLAP evolution in the full kinematic range. 
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Table 1: The measured F2 values. The values of x and at which F2 are determined are 
shown outside the brackets, while the bin boundaries are shown inside the brackets. The number 
of events in each bin as well as the estimated number of background events are given. The 



corrections for Fl, 5l (see eq. |7^, are also given 



x{range) 






F2 i stcit_Jljg 


h% 


= 3.5 


[3.2 - 4.0) GeV" 


6.3(5.0-8.0)10"^ 
1.0(0.8 - 1.3)10"^ 


2196 
3721 


101 
128 


0.875 ± 0.052l[!:°g3 
0.939 ± 0.040l[J:[J?i 


7.8 
2.4 


g2 = 4.5 


[4.0 - 5.0) GeV" 


1.0(0.8 - 1.3)10"^ 
1.6(1.3 - 2.0)10"^ 
2.5(2.0 - 3.2)10"^ 


5131 
4671 
2967 


211 
103 
39 


1.030 ± 0.028±[j:[!i 
1.038 ± 0.0281° 
0.932 ± 0.029l[! °i 


4.5 
1.5 
0.5 


g2 = g,5 (5,0 _ 7.0) GeV" 



1.0 


;o.8 


- i-s; 


10- 


-4 


3192 


179 


1.6 


;i.3 


- 2.0; 


10- 


-4 


8529 


245 


2.5 


;2.o 


-3.2; 


10- 


-4 


9101 


150 


4.0 


;3.2 


-5.0; 


10- 


-4 


6911 


50 


6.3 


;5.o 


-8.0; 


10- 


-4 


5909 


60 


1.0 


;o.8 


- 1.3; 


10" 


-3 


5278 


10 


1.6 


'1.3 


- 2.0^ 


10" 


-3 


4655 


20 



1.221 ±0.036; 
1.138 ±0.021; 
1.038 ±0.018; 
0.951 ±0.019; 
0.841 ±0.018; 
0.776 ±0.017; 
0.753 ±0.019- 



-0.098 
-0.076 
-0.051 
-0.034 
-0.021 
-0.026 
-0.031 
-0.028 
-0.026 
-0.017 
-0.036 
-0.028 
-0.024 
-0.030 



11.8 

3.6 
1.2 
0.4 
0.2 
0.1 
0.0 



= 6.5 (5.0 - 9.0) GeV 



2.5 
4.0 
1.6 
0.4 



2.0 - 3.2 
3.2-8.0 
0.8-3.2 
0.3 - 1.3 



10-3 
10-3 
10-2 
10-1 



9522 

15723 

16255 

1736 



90 
40 
20 




0.645 ± 0.012;0:[]32 
0.629 ± 0.009t[J:°| 
0.495 ± O.OOSttxm 
0.403 ± O.OldtUfs 



0.0 
0.0 
0.0 
0.0 



g2 = 8.5 (7.0 - 9.0) GeV" 



1.6 


1.3 


-2.0) 


10" 


-4 


3666 


109 


1.337 ±0.037: 


hO.060 
-0.081 


7.1 


2.5 


;2.o 


-3.2) 


10" 


-4 


6502 


132 


1.151 ±0.024: 


h0.037 
-0.023 


2.3 


4.0 


;3.2 


-5.0) 


10" 


-4 


6000 


77 


0.993 ±0.021: 


h0.028 
-0.019 


0.8 


6.3 


;5.o 


-8.0) 


10" 


-4 


5983 


40 


0.920 ±0.019: 


h0.014 
-0.025 


0.3 


1.0 


;o.8 


- 1.3) 


10" 


-3 


5827 


50 


0.831 ±0.018: 


h0.031 
-0.025 


0.1 


1.6 


;i.3 


- 2.0) 


10" 


-3 


4657 


30 


0.704 ±0.017: 


h0.023 
-0.015 


0.0 



g2 = 10.0 (9.0 - 11.0) GeV 



1.6 


1.3 


-2.0) 


10" 


-4 


870 


35 


1.359 ±0.075: 


^0.114 
-0.078 


10.9 


2.5 


;2.o 


-3.2) 


10" 


-4 


3904 


106 


1.231 ±0.032; 


f-0.043 
-0.043 


3.4 


4.0 


;3.2 


-5.0) 


10" 


-4 


4415 


57 


1.136 ±0.026: 


f-0.030 
-0.022 


1.1 


6.3 


;5.o 


-8.0) 


10" 


-4 


4479 


28 


0.968 ± 0.022i[':[]30 


0.4 


1.0 


;o.8 


- 1-3) 


10" 


-3 


4370 


48 


0.857 ±0.021: 


f-0.024 
-0.023 


0.1 


1.6 


1.3 


-2.0) 


10" 


-3 


3737 


10 


0.792 ±0.021; 


f0.021 
-0.022 


0.1 


2.5 


;2.o 


-3.2) 


10" 


-3 


3971 


50 


0.715 ±0.018; 


f0.034 
-0.028 


0.0 



Q' = 10.0 (9.0 - 13.0) GeV 



6.3 
1.6 



3.2-8.0 
0.8-3.2 



10-3 
10-2 



12284 
15549 



40 
30 



0.593 ± 0.008l[] []22 
0.517 ±0.007±o;oi3 



0.0 
0.0 
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Table 1 (continued): The measured F2 values. 



x{range) 






7-1 1 ( i-\-S1JS 

r2 ± stat_gyg 


r (7/ 


0.8(0.3-1.3) 


10" 




3685 





0.395 ± 0.010±5:^g§ 


0.0 






Q' 




12.0 (11.0 - 


^ T 

- 13.0) GcV" 




2.5 


;2.0 - 3.2) 


10" 


-4 


2311 


94 


1.300 ± 0.0421^;^^:^ 


5.4 


4.0 


;3.2 - 5.0) 


10" 




3052 


50 


1.128 ± 0.031Io;q24 


1.7 


6.3 


;5.0 - 8.0) 


10" 


A 
-4 


3339 


20 


1.102 ± 0.027t^;^3^ 


0.6 


1.0 


;0.8 - 1.3) 


10" 




3242 


40 


0.903 ± 0.023t[}:E]^9 


0.2 


1.6 


;i.3-2.0) 


10" 


Q 
-0 
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70 


0.789 ± 0.023l^;^^f 


0.1 


2.5 


;2.0 - 3.2) 
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Q 
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10 


0.746 ± 0.021i[];[]^2 


0.0 






Q' 
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- 16.0) GeV^ 




2.5 


;2.0 - 3.2) 


10" 


A 
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84 


1.559 ±0.065lO:?9i 


9.5 


4.0 


;3.2 - 5.0) 


10" 


-4 


3247 


62 


1.338 ± 0.035 


2.9 


6.3 


'5.0 - 8.0) 


10" 


A 


3421 


65 


1.092 ± 0.027ijj;535 


1.0 


1.0 


^0.8 — 1.3^ 


1 n- 
iU 


-3 


3487 


20 


0.992 ± 0.025l[];J];^5 


0.3 


1.6 


1.3 — 2.0^ 


iU 


-3 


2978 


20 


n onfi -4- n no/i +0.021 
U.yUo it U.Uz4_g 


n 1 
U.i 


2.5 


^2.0 — 3.2) 


iU 


-3 


3022 





n 7QQ _i_ fi noi +0.037 
u. 1 c/O in u.uzi_Q Q-|^2 


n n 


6.3 


^3.2 - 8.0) 


1 n- 
iU 


-3 


5535 


10 


n «fio -1- n ni 4+0.020 

U.DD^ It U.Ui4_Q 025 


n n 
u.u 










15.0 (13.0- 


- 20.0) GeV 




1.6 


;0.8 - 3.2) 


10" 




14175 


50 


0.553 ± 0.007l[];^i^ 


0.0 


0.8 


'0.3 - 1.3) 


10" 


1 

1 


5160 





0.370 ± 0.00815;^^^ 


0.0 






Q' 




18.0 (16.0- 


— n 

- 20.0) GeV 




4.0 


;3.2-5.0) 


10" 


-4 


2672 


85 


1.438 ± 0.041j:[};{]^^ 


4.5 


6.3 


;5.0 - 8.0) 


10" 


-4 


3006 


87 


1.166 It 0.031Iq;o37 


1.5 


1.0 


;0.8 - 1.3) 


10" 


Q 
-0 


3118 


39 


1.015 ±0.028l[]:[]i^ 


0.5 


1.6 


;i.3-2.0) 


10" 


Q 
-0 


2682 


29 


0.932 ± 0.0261}]:^^ 


0.2 


2.5 


;2.0 - 3.2) 


10" 


Q 


2737 


10 


0.879 ± 0.024;[];[jff 


0.1 


6.3 


;3.2 - 8.0) 


10" 


Q 



4755 





0.629 ± 0.013i^:^^f 


0.0 






Q' 




22.0 (20.0 - 


^ T^. T 

- 25.0) GcV" 




4.0 


;3.2 - 5.0) 


10" 


/I 


1363 


50 


1.507 ±0.058tg:J^.^ 


7.5 


6.3 


;5.0-8.0) 


10" 


-4 


2568 


60 


1.324 ± 0.037l[J;[]^f 


2.3 


1.0 


;0.8 - 1.3) 


10" 


Q 
-0 


2498 


8 


1.081 ± 0.028Iq;o35[ 


0.8 


1.6 


'1.3 - 2.0) 


10" 


Q 


2136 


20 


0.930 ± 0.0261^:^^1 


0.3 


2.5 


^2.0 — 3.2^ 


lU 


-3 


2240 





0.906 ± 0.0251^26 


0.1 


4.0 


^3.2 — 5.0) 


iU 


-3 


1892 





n 7/in -1- n rm+0.039 
U. /4y it U.Uzz_o 033 


n n 
U.U 


D.O 


'en s n^ 
^O.U — o.U^ 


1 n~ 


-3 


2003 





n 7ni 4- n noi +0.032 

K).iK)l± U.UZi_o 020 


n n 
u.u 






Q 




22.0 (20.0 - 


- 32.0) GeV 




1.0 


;0.8 - 1.3) 


10" 


-2 


3613 





n (^m _L n m q+0.018 
U.OUi it U.Uio_o 030 


U.U 


2.5 


;i.3-5.0) 


10" 


-2 


9701 





0.490 ± 0.007lO;[]i? 


0.0 


0.8 


;0.5 - 1.3) 


10" 


-1 


2325 





0.338 ± 0.009±!):[55^ 


0.0 






Q' 




27.0 (25.0 - 


- 32.0) GeV" 




6.3 


;5.0-8.0) 


10" 


-4 


2314 


60 


1.460 ±0.041tO:[J|^ 


3.8 


1.0 


;0.8 - 1.3) 


10" 


-3 


2356 


35 


1.194 ±0.03210;[]30 


1.2 
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Table 1 (continued): The measured F2 values. 



x{range) 






r2 ± stat_syg 


Ol/o 


1.6 


[1.3 - 2.0; 


10" 


-3 


2060 





1.140 ± 0.033l^^^ 


0.4 


2.5 


^2.0 - 3.2; 


10" 


-3 


1854 


11 


0.928 ± 0.0281^;^^^ 


0.2 


4.0 


;3.2 - 5.0; 


10- 


-3 


1695 


8 


0.765 ± 0.024l^f^ 


0.1 


6.3 


^5.o - 8.0; 


IQ- 


-3 


1709 


2 


0.671 ±0. 02 llJ];J]ji 


0.0 










35.0 (32.0 - 


- 40.0) GeV^ 




6.3 


;5.o - 8.0; 


10-^ 


966 


35 


1.565 ± 0.073l[]:[]^^ 


7.3 


1.0 


^0.8 - 1.3; 


10" 


_3 


1756 


25 


1.370 ± 0.044l[J;[j^| 


2.2 


1.6 


[1.3 - 2.0] 


10" 


_3 


1543 


8 


1.149 ± 0.038Iq;q28 


0.7 


2.5 


[2.0 - 3.2] 


10" 


_3 


1631 





1.018 ± 0.0321^^^ 


0.3 


4.U 


^O.Z — O.U^ 


10" 


_3 


1423 





0.887 ± 0.0301^:5^1 


0.1 


0.0 


o.u — O.U^ 


iU 


-3 


1306 





n 7/1 4- n no7+o.040 
U. ( 4y ± U.UZ ( _Q Qj^3 


n n 
U.U 


1 n 
i.U 


Q 1 
^U.o — i.d^ 




-2 


1307 


10 


f) «oi _)_ n noc+O.OlS 


n n 

U.U 






Q 




35.0 (32.0 - 


- 50.0) GeV 




1.6 


;i.3 - 2.0; 


10" 


_2 


2026 





0.604 ± 0.0171^;^^ 


0.0 


2.5 


;2.o - 5.0; 


10" 


_2 


3845 





0.570 ± 0.012l^;^5^ 


0.0 


0.8 


;o.5 - 1.3; 


10" 




2323 





0.461 ± 0.0131^29 


0.0 






Q' 




45.0 (40.0 - 


- 50.0) GeV^ 




1.0 


'0.8 - 1.3' 


10" 


-3 


1431 


25 


1.441 ± 0.053lg:0^| 


3.9 


1.6 


1.3 - 2.0; 


10" 


-3 


1293 


33 


1.178 ± 0.045Iq;q45 


1.3 


2.5 


^2.0 - 3.2; 


10" 


-3 


1390 


11 


1.071 ± 0.037i^:5^§ 


0.4 


4.0 


;3.2 - 5.0; 


10" 


_3 


1172 


10 


0.894 ± 0.034j:^;5^^ 


0.2 


6.3 


;5.o - 8.0; 


1 A- 

10 


-3 


1129 





0.804 ± 0.031I^;5|| 


0.1 


1.0 


^0.8 - 1.3; 


10" 


_2 


1108 





0.745 ± 0.029lg:{]|| 


0.0 






Q' 




60.0 (50.0 - 


- 65.0) GeV^ 




1.0 


;o.8- 1.3; 


10" 


-3 


790 


8 


1.553 ±0.073l[J:[]^^ 


8.0 


1.6 


;i.3-2.o; 


10" 


-3 


1250 


8 


1.339 ± 0.050lE]:[]il 


2.5 


2.5 


;2.o - 3.2; 


10" 


-3 


1294 


8 


1.100 ±o.o4oi[;:[]i 


0.8 


4.0 


^3.2 - 5.0^ 


10" 


-3 


1166 





1.012 ± 0.039l^:^4f 


0.3 


6.3 


;5.o - 8.0; 


10" 


Q 
-0 


1118 





r\ or\r\ i r\ Aor-l-n 044- 

0.890 ± 0.035Io;q2o 


r\ -1 

0.1 


1.0 


'0.8 - 1.3' 


10" 


-2 


952 





0.672 ± 0.028lg:g?| 


0.0 


1.6 


;i.3 - 2.0; 


10" 


-2 


760 





0.614 ± 0.029lE]:[J|| 


0.0 


2.5 


;2.o - 5.0; 


10" 


-2 


1719 





0.627 ± 0.020j:[]:{]?[] 


0.0 


0.8 


;o.5 - 1.3; 


10" 


-1 


1391 





0.521 ± 0.0191^:^1^ 


0.0 






Q' 




70.0 (65.0 - 


- 85.0) GeV" 




1.6 


;i.3 - 2.0; 


10" 


-3 


862 


25 


1.318 ± 0.062l[]:[]66 


3.6 


2.5 


^2.0-3.2; 


10" 


-3 


1107 


20 


1.228 ± 0.0491[]:{]^^ 


1.2 


4.0 


;3.2 - 5.0; 


10" 


-3 


885 





0.917 ±0.039l[J:[]f3 


0.4 


6.3 


;5.o-8.o; 


10" 


-3 


852 


10 


0.887 ± 0.040l[!:[]i? 


0.1 


1.0 


;o.8- 1.3; 


10" 


-2 


816 





0.774 ± 0.035l[]:0i^ 


0.0 


1.6 


;i.3-2.o; 


10" 


-2 


616 





0.697 ± 0.0371[]:[]36 


0.0 


2.5 


;2.o-5.o; 


10" 


-2 


1390 





0.607 ±0.0211[]:0^9 


0.0 
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Table 1 (continued): The measured F2 values. 



x{ra,nge) 


A T 

Nev 


A T 


P2 ± stat^sys 


r (7/ 


0.8(0.5-1.3; 


10" 


_\ 


1244 





0.484 ± 0.018io:o26 


0.0 










90.0 ( 


35.0 - 


110.0) GeV^ 




1.6 


;i.3-2.o; 


10- 




299 





1.363 iO.lOllO;!!^ 


6.6 


2.5 


;2.o-3.2; 


10" 




746 


8 


1.275 ± 0.063l[]:[]i 


2.1 


4.0 


;3.2 - 5.0; 


10" 


-3 


702 


8 


1.068 ± 0.0541[]:046 


0.7 


6.3 


;5.o-8.o; 


10" 


Q 

-0 


689 





0.974 ± 0.049lE]:[ji 


0.2 


1.0 


;o.8- 1.3; 


10" 


-2 


667 





r\ r-\ A ^ 1 r\ r\ A i-\_l_n 0^*3 

0.845 ± 0.043l[J:[]p 


0.1 


1.6 


;i.3-2.o; 


10" 


-2 


456 





0.619 ± 0.037lo:o58 


0.0 


2.5 


;2.o - 3.2; 


10" 


-2 


468 





0.556 ± 0.033l5:^fg 


0.0 


4.0 


;3.2 - 5.0; 


10" 


-2 


473 





0.526 ± 0.0311[]:[]|^ 


0.0 


0.8 


;o.5- 1.3; 


10- 


-1 


987 





0.453 ± 0.019l[j;[3f^ 


0.0 


2.0 


;i.3-3.2; 


10- 


-1 


122 





0.255 ± 0.029lg;^§^ 


0.0 




Q'- 


= 120.0 (110.0 


- 140.0) GeV^ 




2.5 


;2.o - 3.2; 


10- 


-3 


477 


8 


1.467 ±0.09310:038 


4.1 


4.0 


;3.2-5.o; 


10- 


-3 


519 


17 


1.176 ±0.073l[J-030 


1.3 


6.3 


;5.o-8.o; 


10- 


-3 


520 





1.112 ±0.066lo;oi 


0.4 


1.0 


;o.8- 1.3; 


10- 


-2 


466 





0.836 ± 0.0511^;[]f^ 


0.1 


1.6 


;i.3-2.o; 


10- 


-2 


321 





0.610 ± 0.043Iq;q23 


0.0 


2.5 


;2.o - 3.2; 


10- 


-2 


366 





0.653 ± 0.0451^;^^^ 


0.0 


4.0 


;3.2- 5.0; 


10- 


-2 


345 





0.519 ± 0.036T^;^if 


0.0 


0.8 


;o.5- 1.3; 


10- 


-1 


716 





0.451 ± 0.022l^;^ig 


0.0 


2.0 


;i.3-3.2; 


10- 


-1 


125 





0.207 ± 0.022±g;[]|f 


0.0 




Q'- 


= 150.0 (140.0 


- 185.0) GeV^ 




2.5 


^2.0-3.2; 


10- 


-3 


159 





1.275 ±0.130lE]:[jp 


7.1 


4.0 


;3.2-5.o; 


10- 


-3 


376 


8 


1.013 ± 0.07010:026 


2.1 


6.3 


;5.o-8.o; 


10- 


-3 


418 





0.991 ± 0.063l[]:[]i 


0.7 


1.0 


;o.8- 1.3; 


10- 


-2 


382 





0.835 ± 0.056l[]:[]?6 


0.2 


1.6 


;i.3-2.o; 


10- 


-2 


316 





0.685 ± 0.050lo;o^3 


0.1 


2.5 


^2.0-3.2; 


10- 


-2 


314 





0.719 ± 0.0541^;}]^^ 


0.0 


4.0 


;3.2-5.o; 


10- 


-2 


297 





0.522 ± 0.039l[]:[]}8 


0.0 


0.8 


;o.5- 1.3; 


10- 


-1 


573 





0.420 ± 0.022l[J;[jJ^ 


0.0 


2.0 


;i.3-3.2; 


10- 


-1 


204 





0.288 ± 0.026l[]:24i 


0.0 




Q'- 


= 200.0 (185.0 


- 240.0) GeV^ 




4.0 


;3.2-5.o; 


10- 


-3 


186 





1.035 ± 0.09610:026 


4.2 


6.3 


;5.o-8.o; 


10- 


-3 


267 





1.041 ± 0.08410-027 


1.3 


1.0 


;o.8- 1.3; 


10- 


-2 


258 





0.832 ± 0.067lg:g|^ 


0.4 


1.6 


1.3-2.0; 


10- 


-2 


230 





0.720 ± 0.06210-021 


0.1 


2.5 


;2.o - 3.2; 


10- 


-2 


189 





0.631 ± 0.06010-032 


0.0 


4.0 


;3.2- 5.0; 


10- 


-2 


194 





0.541 ± 0.05010-0^3 


0.0 


0.8 


;o.5- 1.3; 


10- 


-1 


372 





0.410 ± 0.02710-016 


0.0 


2.0 


;i.3-3.2; 


10- 


-1 


178 





0.271 ± 0.02510:012 


0.0 
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Table 1 (continued): The measured F2 values. 



x{range) 






7-1 1 ( j.-\-S1JS 

F2 ± statl/y^ 


Sl% 










= 250.0 (240.0 


^ T 

- 310.0) GeV^ 




4.0 


;3.2 


— 


5.0)10" 


-3 


67 





1.353 ± 0.219t[!:[3|4 


7.2 


6.3 
1.0 


;5.o 
;o.8 


— 


8.0)10" 
1.3)10" 


-3 
-2 


184 
200 






1.176 ±0.11610:[]30 
0.916 ± 0.0851^:^26 


2.1 
0.6 


1.6 


;i.3 


— 


2.0)10" 


-2 


147 





0.699 ± 0.074Io;q^3 


0.2 


2.5 


;2.o 


— 


3.2)10" 


-2 


127 





0.544 ± 0.0611E]:[]i^ 


0.1 


4.0 


;3.2 


— 


5.0)10" 


-2 


148 





0.583 ± 0.063j:[J;[3^^ 


0.0 


0.8 


;o.5 




1.3)10" 


-1 


275 





0.435 ± 0.034j:[];[}}g 


0.0 


2.0 


;i.3 


— 


3.2)10" 


-1 


148 





0.267 ± 0.0271[]:[]}J 


0.0 










= 350.0 (310.0 


- 410.0) GeV^ 




6.3 
1.0 


;5.o 
;o.8 


— 
— 


8.0)10" 

1.3)10" 


-3 
-2 


129 

144 







1.107 ±0.127±g:0|4 

0.814 ± 0.086t[];[]f^ 


4.8 

1.3 


1.6 


1.3 


— 


2.0)10" 


-2 


120 





0.661 ± 0.076Iq;q];9 


0.4 


2.5 


;2.o 


— 


3.2)10" 


-2 


95 





0.446 ± 0.055To;oo9 


0.1 


4.0 
0.8 


;3.2 
;o.5 


— 
— 


5.0)10" 
1.3)10" 


-2 
-1 


103 
199 






0.672 ± 0.089l^;^^g 
0.414 ± 0.0371^}'^ 


0.0 
0.0 


2.0 


;i.3 


— 


3.2)10" 


-1 


130 





0.262 ± 0.029l{];[J(^,^^ 


0.0 








Q'- 


= 450.0 (410.0 


- 530.0) GeV^ 




6.3 


;5.o 


— 


8.0)10" 


-3 


51 





1.722 ± 0.34810:[]82 


8.7 


1.0 


;o.8 


— 


1.3)10" 


-2 


89 





0.929 ± 0.128l[!:[]|[j 


2.4 


1.6 


;i.3 


— 


2.0)10" 


-2 


84 





0.864 ± 0.125Iq;q33 


r\ '~7 

0.7 


2.5 


;2.o 




3.2)10" 


-2 


89 





0.688 ± 0.095lo:o26 


0.2 


4.0 


;3.2 


— 


5.0)10" 


-2 


52 





A c\ 1 n\ /'A r-r _l_0 01 7 

0.435 ± 0.075l^:5}4 


0.1 


0.8 


;o.5 


— 


1.3)10" 


-1 


148 





0.455 ± 0.048l[]:[]i? 


0.0 


2.0 


;i.3 


— 


3.2)10" 


-1 


115 





0.348 ± 0.0431^;^}^ 


0.0 








Q'' = 


= 650.0 (530.0 


^ !T 

- 710.0) GeV^ 




1.0 


;o.8 


— 


1.3)10" 


-2 


76 





1.067 ± 0.163l[!:g^| 


5.7 


1.6 


1.3 


— 


2.0)10" 


-2 


71 





0.897 ± 0.142Iq;q29 


1.6 


2.5 
4.0 


;2.o 
;3.2 


— 
— 


3.2)10" 
5.0)10" 


-2 
-2 


78 
43 




10 


0.728 ± 0.1091^;^^^ 
0.471 ± 0.102l^;^}| 


0.5 
0.1 


0.8 


;o.5 


— 


1.3)10" 


-1 


104 





0.445 ± 0.0561^:^20 


0.0 


2.0 


;i.3 


— 


3.2)10" 


-1 


96 





0.344 ± 0.046l^:^}4 


0.0 








Q'- 


= 8 


00.0 (710.0 


- 900.0) GeV" 




1.0 


;o.8 




1.3)10" 


-2 


22 





0.890 ± 0.2411[]:[j|| 


9.3 


1.6 


;i.3 




2.0)10" 


-2 


40 





0.819 ±0.169l[!-[]29 


2.6 


2.5 


;2.o 




3.2)10" 


-2 


48 





0.657 ± 0.122 w 


0.8 


4.0 


;3.2 




5.0)10" 


-2 


32 





0.512 ±0.115l[;:[ji9 


0.2 


0.8 


;o.5 




1.3)10" 


-1 


54 





0.293 ± OMSt^fol 


0.0 


2.0 


;i.3 




3.2)10" 


-1 


45 





0.227 ± 0.04110:0$8 


0.0 
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Table 2: The measured F2 and F|™ values for Q'^ > 1000 GeV^ (see text). The values of 
X and at which F2 and F^^"* are determined are shown outside the brackets, while the bin 
boundaries are shown inside the brackets. The number of events in each bin as well as the 
estimated number of background events are given. The corrections for F^, 5l (see eq. [7^, are 
also given. 



x{range) 






F2 


F2 i sta,t_gyg 








Q' 




1200.0 (900.0 - 1300.0) GeV^ 




2.5(2.0- 


2.0 
3.2 


)io- 
)io- 


-2 
-2 


34 
41 






0.645 
0.652 


0.634 ±0.135;0:[]30 
0.641 ±0.1281° 


6.8 
1 9 


4.0(3.2- 


5.0 


)io- 


-2 


32 





0.586 


0.577 ±0.132l[J 


0.5 


0.8(0.5- 
2.0(1.3- 
5.1(3.2- 


1.3 
3.2 
8.0 


)io- 
)io- 
)io- 


-1 
-1 
-1 


54 
48 
11 







0.321 
0.308 
0.056 


0.316 ± 0.052j:0:[]J2 
0.304 ± 0.057;[] 
0.055 ± 0.020t[j^[}i 


0.1 
0.0 
0.0 










1500.0 (1300.0 - 1800.0) GeV^ 




2.5(2.0- 


5.0 


)io- 


-2 


52 





0.614 


0.601 ±0.105t[];[]i 


3.2 


0.8(0.5- 


1.3 


)io- 


-1 


46 





0.605 


0.593 ±0.119t[]:[ji 


0.1 
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Figure 1: Distance (in cm) between track and RCAL position a) in X, b) and Y. Distance (in 

cm) between track and SRTD position c) in X, d) and Y. e) Vertex reconstruction efficiency as 
a function of (points: data; histogram: MC). f) Vertex resolution along the beam (in cm) 
as a function of 7^ . 
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Figure 2: The ratio averaged over prh a* a function of for several ranges of'y^. The 

curves show the correction function C{^^,pTh, In) > o-veraged over pth a function of for 
several ranges in 7^ . 



30 




Figure 3: a) The distributions of (dashed- dotted), (dashed), (dotted) and 

(solid) for several ranges in 7^ as determined from the MC simulation, b) The distri- 
butions of (dotted), (dashed) and (solid) for several ranges in 7^ as obtained 

from the MC simulation, c) The corresponding distributions of (solid), §^ (dashed), 
(dotted) and (dashed- dotted). 
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Figure 4: The distributions of ^ for a) 140° < 7^ < 160° and h) 160° < 7^ < 180°. The 
points show the data and the shaded histograms show the results from the MC simulation. 
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Figure 5: Distributions of a) positron energy, b) positron angle, c) Zyertex, d)-f) ^ for different 
ranges of'y^. The vertical lines in c)-f) indicate the positions of the cuts used in the analysis. 
For d)-f) the cut values are shown for the central value of 'y^ in each bin (see text). The MC 
distributions are normalised to the integrated luminosity of the data. 
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Figure 6: The distributions of a) yp^, h) Xp^, c) Q^p^and d) 7^^, for data (points) compared with 
the distributions obtained from the MC simulation (shaded histograms). The MC distributions 
are normalised to the integrated luminosity of the data. 
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Figure 7: a)The distribution of the events in the (x,Q'^) plane, b) The (x,Q'^)-bins used in the 
structure function determination. Also indicated are lines of constant y and of constant 7^. 
The 7^ values for this figure are calculated directly from x and . 
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Figure 9: The distribution of Tn{Ei — pzi) for data (solid dots), DIS MC simulation (solid 
histogram), photoproduction (PHP) MC simulation (dashed histogram) and LUMI tagged pho- 
toproduction data (open squares). The PHP distribution is valid for T,i{Ei — pzi) > 38 GeV. 
The result of the fit to the data, described in the text, is shown as the dashed-dotted curve and 
the photoproduction contribution to this fit is indicated by the dotted curve. 
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Figure 10: Relative systematic errors as a function of y for different categories of systematic 
uncertainties for bins with < 100 GeV^. The total systematic and statistical errors are 
also given. The bottom two plots show total systematic and statistical errors for bins with with 
> 100 GeV^. In categories 4) (ind 5), and the total systematic error plot for < 100 
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Figure 11a: Structure function F2 as a function of x for fixed values of Q"^ (1.5 GcV^ < < 
22 GeY'^J.The solid dots correspond to the data from this analysis. The "ISR" (open squares) 
and "SVX" (open triangles) data show the ZEUS results obtained previously. The inner thick 
error bars represent the statistical error, the full error bars correspond to the statistical and 
systematic errors added in quadrature. The results from NMC are also shown (open circles). 
The solid lines indicate the QCD NLO fit to the data used for the acceptance correction. The 
values are indicated in units of GeV^ . 
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Figure 11c: Structure function F2 (F^"" for > 1000 GeV^) as a function of x for fixed values 
ofQ^ (450 GeV^ < < 5000 GeV^). The values are indicated in units of GeV^. 



41 



ZEUS 94 



1:^ 



• Thi;s analysis 

// / □ ZEUS ISR 
ilj A ZEUS SVX 

-/ ^ X= 6.3E-05 




— ZEUS EIT 

I'' MRSA' 

/ ; y GRV94 
/ ;]i CTEQ 3 

' U JIC= 0.0001 

4 




" // X- 0.0001 6 




: # 

- ,jK Jic= 0.00025 




^ X= 0.0004 i 


X= O.OG06.5 , 


X- 0.00'; 


X= 0.001 6 










- M 

1 1 1 1 1 1 1 1 


''''^^ 1 1 1 







10 



10 



10 



- JIC= 0.0025 

1 1 1 


- X= 0.004 

1 1 1 


- JC= 0.0063 

1 1 1 1 


- X= 0.01 

1 1 1 1 


\ X= 0.01 6 

1 1 


\ X= 0.025 

1 1 1 1 


JC= 0.04 

1 1 1 1 


~ X= 0.081 


- X= 0.2 

1 1 


- JC= 0.51 

1 1 1 1 1 !• r~ 


10 10^ 10^ 10 10^ 10^ 



10 10^ lO'^ 



10 10^ lO'' 



Figure 12: Structure function Fl'"* as a function of for fixed values of x. The data from this 
analysis correspond to the solid dots. The "ISR" (open squares) and "SVX" (open triangles) 
data show the ZEUS results obtained previously. The solid line indicates the QCD NLO fit to 
the data used for acceptance correction. Also indicated are the MRSA', GRV94 and CTEQ3 
parameterisations. 
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Figure 13a: Structure function F2 as a function of x for fixed values of . The solid dots 
correspond to the data from this analysis. The "ISR" (open squares) and "SVX" (open trian- 
gles) data show the ZEUS results obtained previously. The inner thick error bars represent the 
statistical error, the full error bars correspond to the statistical and systematic errors added in 
quadrature. The results from NMC, E665, BCDMS and SLAC are also shown (open circles). 
The solid line indicates the QCD NLO fit to the data used for acceptance correction. Also 
indicated are the MRSA' , GRV94 and CTEQ3 parameterisations. The values are given in 
GeV2. 
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Figure 13b: Structure function F2 as a function of x for fixed values of Q^. The solid line 
indicates the QCD NLO fit to the data used for acceptance correction. Also indicated are the 
MRSA' , GRV94 and CTEQ3 parameterisations. The values are given in GeV^. 
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A Systematic uncertainties and their correlations 



In Sect. |7]1|, the determination of the systematic uncertainties of the F2 values was discussed. 
The resulting systematic errors are given in Tables |l] and 0. Bin to bin correlations exist for 
the systematic errors. Given the small statistical and systematic uncertainties of the present 
measurement these correlations should be taken into account when performing fits. 



A.l F2 values resulting from systematic checks 



Tables |Ala| and |Alb| give the F2 values obtained for each of the 29 different checks of systematic 



uncertainties described in Sect. ^Tll fl As noted above, the systematic uncertainties for the F2 
measurements at > 100 GeV^ and 1/ = ~" > 0.01 are small compared to the statistical 
errors. Therefore, these F2 values are given only for the data with < 100 GeV^ or y < 0.01. 

In the calculation of the total systematic uncertainty given in Tables |l] and 0, the deviations 
from the central value for each check were added in quadrature, except for the two checks 
concerning diffractively produced events (n" 23, 24) whose contributions were calculated from 
the maximum offset, positive or negative, produced by the checks of n" 23 and 24. 

Fits to the F2 data, which take the systematic error correlations into account, may be performed 
in the following way. First the F2 values with their statistical errors as given in Tables ^ and |] 
are fitted. This yields the standard values of the fit parameters (oj) and their statistical errors 
{ai{stat)). The fit is then repeated for each of the 29 sets of F2 values listed in Tables |Ala 



and |Alb| resulting in 29 different sets of fit parameters af , / = 1 — 29. An estimate of the 
systematic errors of the fit parameters may be obtained by adding in quadrature, for each fit 
parameter, the differences from the central parameter values, ai{syst) = \Jj2fio,{ ~ (^lY- 

Note that the overall normalisation error of 2% (see Sect. [7.1|) must be taken into account 
separately. 



A. 2 Compact parametrisation of uncertainty correlations 



Although we recommend the use of Tables [Ala] and |Alb| , in determining the effect of correlated 
uncertainties in a fit to F2, the large number of systematic checks involved may make this a 
cumbersome process. Therefore, a somewhat more compact representation of the correlated 
systematic uncertainties is provided below. 

We have investigated the dependence on y and of the relative deviations, A{F2)/F2, of the 
F2 values listed in Tables |Ala| and [Alb| from the central values of F2. These deviations can be 



^ This Table and a short description of its recommended usage is available in electronic form from the ZEUS 

It can also be obtained by contacting the 



WWW page whose current address is http://zow00.desy.de:8000/ 
authors. 
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grouped in six different classes of systematic behaviour, one of which is a function of only 
and the others are functions of y only. 



A ■ 
B : 
C : 
D : 
E : 
F : 



A(F2 



F2 
A{F2 



F2 
A{F2 



F2 
A(F2 



F2 
A(F2 



F2 
A(F2 



F2 



aA-bA- logio 
(^b/\ Iogio2/P 

ac/\ logioZ/P -bc-\ logio?/!^ 
^d/I logiol/P + bD-\ logiol/r - CD 
as ■ I logiol/r - ■ I logio2/|^ 
-ap + 6f ■ I logioZ/l 



The deviations produced by the systematic checks fall into one of the classes listed above 
(positive class) or into a class (negative class) where the signs of the coefficients are reversed 
(these classes are denoted by, e.g. -A). The parameters are determined by fitting the deviations 
produced by each check separately to the appropriate functional form. For each class the 
parameters obtained from each fit are then added in quadrature. 



The results for the six classes of systematic errors are given in Table ^ and their dependence 
on and y are displayed graphically in Fig. Note that these functions only represent 

the systematic errors within the kinematic range of the present F2 measurements. 

These functional forms may be used for systematic studies to reproduce the offsets of the values 
of F2 for < 100 GeV^ or y < 0.01 given in this paper. Outside this kinematic region, the 
systematic errors given in Table |T] can be treated as uncorrelated errors. 

The systematic deviations of the F2 values as described by the 6+6 classes determined can be 
added to the central values of F2 to produce 12 sets of F2 data. These 12 sets may be used. 



in the same way as described in Sect. |A.1| for the 29 sets, to determine the effect of correlated 
uncertainties on fits to F2. 

These functions have been tested using the QCD NLO fit described in Sect. ^ The systematic 
uncertainties of the fit parameters are consistent with the more complete treatment described 
in Sect. |A.1| . It should be noted, however, that the simplification of the systematic correlations 
results in some loss of information, and that the use of the full treatment described in Sect. [A.l| is 
recommended if the most reliable estimates of the effect of correlated uncertainties are required. 



^ A FORTRAN program of these functions is available through the ZEUS WWW home page whose current 
address is http://zow00.desy.de:8000/. It can also be obtained by contacting the authors. 
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Table Ala: The extracted F2 values for each systematic check at values of < 100 GeV^ or y = — < 0.01. 
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Table Alb: The extracted F2 values for each systematic check at values of < 100 GeV^ or y — — < 0.01 



Class 


Checks contri- 
buting to class 


^class 


bclass 


C-dass 


A 


11;5;7 


2.34 


0.96 





-A 


6;8;9 


-1.94 


-0.83 





B 


1;29 


0.126 








B 


2;28 


-0.085 








C 


12;14 


0.153 


4.4-10-3 





-C 


3;4;13;15;27 


-0.145 


-3.4-10--'' 





D 


10;16;26;25 


0.132 


4.3-10-3 


0.69 


D 


17;22;24;23 


-0.146 


-7.4-10-3 


-1.29 


E 


21 


1.62 


0.39 





-E 


20 


-3.05 


-0.77 





F 


19 


1.05 


0.83 





-F 


18 


-1.00 


-0.79 






Table A2: The parameters of the functions describing the different classes of systematic un- 
certainties of F2 (see text). 
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Figure Al: The or y dependence of for the six classes of systematic uncertainties (see 
text). 
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